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Abstract 



We report on recent results concerning the derivation of effective evolution equa- 
tions starting from many body quantum dynamics. In particular, we obtain rigorous 
derivations of nonlinear Hartree equations in the bosonic mean field limit, with precise 
bounds on the rate of convergence. Moreover, we present a central limit theorem for the 
i " fluctuations around the Hartree dynamics. 

^ ; 1 Introduction 

We consider systems of N particles, which can be described in quantum mechanics by a 
wave function ip^ £ L 2 {E? N , dx± . . . dx^)- Here x\ . . . , x^ £ M 3 parametrize the positions 
of the N particles. In nature, one distinguish between bosonic and fermionic systems, whose 
wave functions are characterized by different permutation symmetries. In these notes, we 
will focus exclusively on bosonic systems; correspondingly, we will always assume tpN to be 
symmetric with respect to permutations of the N particles. 



The wave function ipN has a probabilistic interpretation. More precisely, the absolute 
value \iPn(xi, . . . ,xn)\ 2 is interpreted as the probability density for finding the N particles 
CN ■ close to x\, . . . , xn ■ According to this probabilistic interpretation, we will always assume the 

normalization H^jvlh = 1- Apart describing the distribution of the iV particles in space, the 
wave function ipN also determines the probability law of every other observable. A physical 
observable in quantum mechanics is a self-adjoint operator B on L 2 (M. 3N ). If, for example, 
B = J2j \7'l 7r j)( 7r jl> f° r a sequence of real numbers and an orthonormal basis {iTj}j of 

L 2 (M 3JV ), the expectation of B in the state described by the wave function ipx is given by 
the L 2 -product 

(ip N ,B4; N ) = \(xj,i/)N) I 2 • 

3 

This formula is interpreted as follows; the physical observable B assumes the value \j with 
probability iPn)\ 2 - The normalization of the wave function guarantees that 



Y^\(^^n)\ 2 = Un\\1 = 1. 
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A similar interpretation can be easily obtained for observables with continuous spectrum. 
The fact that the same wave function determines the probability law of all physical observ- 
ables implies that measurements of non-commuting observables are not independent. This 
is basis for Heisenberg's uncertainty principle. 

The time evolution of an N particle wave function ip^ G L 2 (M. 3N ) is governed by the 
many-body Schrbdinger equation 

id t ipN,t = H N ip N j ■ (1) 

The subscript t indicates the time-dependence of the wave function. On the right hand side 
of (HJ), Hn is a self-adjoint operator on the Hilbert space L 2 (M. 3N ), known as the Hamilton 
operator of the system. We will consider Hamilton operators of the form 

N N 
j=l i<j 

where V models the interactions, and A £ 1 is a coupling constant (one could also intro- 
duce an external potential; the analysis presented in these notes would still apply). The 
Schrbdinger equation ([1]) is a linear partial differential equation. It always has a unique 
global solution, which can be obtained by applying the unitary group generated by Hn to 
the initial wave function ipN,t=o, i-e. ipN,t = e~ tHNt ipN,o- Hence, the well-posedness of the 
many-body Schrbdinger equation is never an issue. What makes the study of the Schodinger 
equation challenging from the mathematical point of view is the fact that, in typical systems 
of interest in physics and chemistry, the number of particles N involved in the dynamics 
is very large, varying from N ~ 10 3 for extremely dilute Bose Einstein condensates, up to 
values of the order N ~ 10 23 for typical samples in chemistry. For these huge values of N, 
it is in general impossible to extract useful qualitative or quantitative information from (TjQ), 
going beyond the mere existence and uniqueness of the solution. For this reason, one of 
the main goals of non-equilibrium statistical mechanics is the derivation of simpler effective 
evolution equation which provide a good approximation of (pQ) in the interesting regime. 

A simple regime where such effective equations can be derived is the so called mean field 
limit of quantum mechanics. The mean field regime is characterized by a large number of 
very weak collisions among the particles. This limit can be realized by choosing 

N > 1, A < 1, with NX oil fixed. 

The condition N 3> 1 guarantees that each particle experiences a large number of collisions; 
A <C 1 implies, on the other hand, that each collision is very weak. Finally, NX of order 
one means that the total force produced by the many weak collisions is of order one, and 
therefore comparable with the inertia of the particle. To study the time evolution in the 
mean field regime, we can therefore analyze the dynamics generated by the Hamiltonian 

N N 

H N = J2 -A,, + - J2 V( Xl - Xj ) (2) 

j = l i<j 

in the limit of large N. In particular, we are interested in the evolution of initially factorized 
states. We assume therefore that, at time t = 0, 

N 

^JV,0 = ¥® N (meaning that tpN,o(xi, x N ) = <p(xj)) 

3=1 
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for some <p G L 2 (M 3 ). Because of the interaction, factorization is not preserved by the 
time-evolution. Still, considering the mean-field nature of the interaction, we may expect 
factorization of the evolved many body wave function to be approximately restored in the 
limit of large N. In other words, we may expect that, in an appropriate sense, 

1>N,t * <P? N (3) 

as N — > oo. If we assume that this is indeed the case, it is simple to derive a self-consistent 
equation for the evolution of (ft- I n fact, (approximate) factorization of ipN,t implies that 
the iV particles are distributed independently in space, with the common probability density 
function |<^t| 2 . This means that the total potential experienced, say, by particle j, can be 
approximated by 




Hence, if factorization is preserved in the limit of large N, (ft must evolve according to the 
self-consistent nonlinear Hartree equation 

idm = -Atpt + (V* \<pt\ 2 )(ft (4) 
where the many-body interactions have been replaced by an average, mean field, potential 

{v*Wt\ 2 )- 

In which sense can we expect ([3]) to hold true? To answer this question, we introduce 
the reduced density matrices associated with the solution of the Schrodinger equation tpN,t- 
For k = 1,2, . . . , N , we define the /c-particle reduced density as 

7N,t = tv k+l,...,N \ipN,t){^N,t\ 
where trfc+i,...^ denotes the partial trace over the last (N—k) particles, and where |V ; iV,t)(V ; -/v,t| 

(k) 

denotes the orthogonal projection onto ViV.t- In other words, ^y Nt is defined as the non- 
negative trace class operator on L 2 (IR 3fc ), with the kernel 

i k ) ( i ' \ 

7N,t\ x ^ > • • • ' x ki x li ■ ■ ■ i x k) 

= J dx/c + i . ..dx N ipN,t(xi, ■ ■ .,x k ,x k+1 , . . . ,x N )ip N j(x[, . . .,x' k ,x k+1 , . . . ,x N ). 

(k) 

From the normalization of if)N,t, we find tr j Nt = 1 for all k,t,N. Of course, for k < N, 

(k) 

the fc-particle density 7^ does not contain the full information about the state described 

(k) 

by tpN,t- Nevertheless, ^ Nt is sufficient to compute the expectation of an arbitrary A;-particle 
observable. In other words, 

{^N,t, (J W O l (N - k) ) *N,t) = tr (J (fe) (8 l (iV - fc) ) l^tX^tl = tr 7 W 

for any self-adjoint operator JW on L 2 (M ik ). It turns out that reduced density matrices 
provide the right language to describe the convergence of the full many-body Schrodinger 
evolution generated by ([2]) towards the limiting mean field Hartree dynamics Q. Under 
appropriate assumptions on the potentials, we have the following theorem. 
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Theorem 1.1. Assume that, at time t = 0, ipN = , for ip G H 1 (M 3 ). Let ij)N,t = 
e~ tHNt ipN be the solution of the Schrddinger equation with mean field Hamiltonian (dj) ; and 
with initial data ipN- Then, for every fixed k £ N, t £ R, we have 

as N — > oo, with respect to the trace class topology. Here ipt is the solution of the nonlinear 
Hartree equation with initial data (pt=o = f- 

This theorem explains in which sense ([3]) should be understood; the reduced densities of 
ipN,t converge, as N — > oo, towards the reduced densities of the product ff N - In terms of 
expectation of observables, the theorem implies that, for any bounded self-adjoint operator 
on L 2 (R 3k ), 

as N — > oo. In other words, if we compute the expectation of an observable depending only 
on a fixed number of particles k then, in the limit of large N, we can replace the solution of 
the full many-body Schrodinger equation with products of the solution of the one-particle 
Hartree equation. 

The first proof of Theorem 11.11 was obtained in [8] for the case of bounded potentials 
|| V||oo < oo- The techniques introduced in [8] were then extended in [3] to potentials with a 
Coulomb singularity V(x) = ±l/|x|. In [7], a different approach was then developed to show 
the convergence towards the Hartree dynamics for potentials with Coulomb singularities. 
In contrast with the previous methods, this new approach provided precise bounds on the 
rate of the convergence towards the limiting Hartree dynamics. More recently, a different 
derivation, covering also potentials with more severe singularities, was introduced in [B]. A 
more complete list of related works can be found in [21 [1]. In the rest of the paper, the 
approach developed in [7] will be presented, and some of its applications will be discussed. 

2 Fock space representation and evolution of coherent states 

We define the bosonic Fock space over L 2 (M 3 ) as 

T = C © L 2 s (R 3n , 6x x ... dx n ) , 

n>l 

where Lg(M 3n ) denotes the subspace of L 2 (M 3n ) consisting of permutation symmetric func- 
tions (in accordance with the bosonic symmetry). Elements of J- are sequences ip = 
n >i, where ip( n ^ G L 2 (M. 3n ) is an n-boson wave function. J 7 is a Hilbert space with 
respect to the inner product 

^;^)=^ (0) + E^ (n) ^ (n) )- 

n>l 

The idea beyond the construction of the Fock space is that, on J 7 , one can describe states 
with variable number of particles. The state described by the sequence n >o has 

n particles with probability HV^Hl (the normalization = 1 is always imposed). 

To measure the number of particles in a state tp, we introduce the number of particles 
operator M defined by 

(AAV) (n) = nV> (n) ■ 
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Eigenvectors of M are sequences of the form {0, 0, . . . , 0, tp( n \ 0, . . . } having only one non- 
zero component. A special example of such a state is the vacuum f2 = {1,0,...}, describing 
a state with no particles at all. 

In order to define a time-evolution on J 7 , we introduce the Hamilton operator Hn, defined 

by 



N 

n i n 

j=l i<j 



(5) 



By definition, the Hamilton operator T~Ln leaves all sectors with fixed number of particles 
invariant. In particular, when applied on states with exactly N particles, the Hamiltonian 
%M acts exactly as the mean field Hamiltonian ([2]). 

On T it is useful to introduce creation and annihilation operators. For / 6 L 2 (IR 3 ), we 
define namely the creation operator a*(f) and the annihilation operator a{f) by the formulas 

1 n 

(a*(/)V^) (n) (xi, • • • , x n ) = -= f{xM {n ' 1] ^U 'V <'„) 

(a(/)V;) (n) (xi,...,x„) = Vn + l / dxJi^^ n+1 \x, Xl ,...,x n ). 



It is simple to check that creation and annihilation operators are adjoint to each other. Also, 
they satisfy the canonical commutation relations 

[a(f),a*(g)] = (f,g) L , Hf),a(g)] = [a* (/) , a* (g)] =0. (6) 

Applying creation operators on the vacuum, it is possible to generate an orthonormal basis 
for the Hilbert space T . For example, the TV-particle state tp® N can be expressed as 

{0,...,0,^ iV ,0,...} = ^L=(a*M)^. 
We will also make use of operator valued distributions a*,a x , defined by 
«*(/) = J dxf(x)a* x a(f) = J dx f(x) a x . 

In terms of these operator valued distributions, we can express the number of particles 
operator as 

N = I dx a* x a x . 



A simple consequence of this formula is that, although creation and annihilation operators 
are not bounded, they are nevertheless bounded with respect to the square root of the number 
of particle operator, in the sense that 

IK/MI < imhHAA 1 / 2 ^ and \\a*Um<\\fhW+l? ,2 n (7) 
for all / G L 2 (R 3 ), ip £ T. 



Also the Hamilton operator T~Ln can be expressed in terms of the operator valued distri- 
butions 

V.N = J dxV xC&V x a x + J dxdyV (x - y)a x a*a y a x . (8) 
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The fact that the Hamiltonian commutes with N is reflected, in (|8|), by the observation that 
each term in T~Ln contains the same number of creation and annihilation operators. As a 
consequence the evolution generated by Hn leaves the sectors of T with fixed number of 
particles invariant. In particular, if we consider the evolution of an initially factorized state, 
we find 

e- iHNt {0, . . . , 0, v? 07V , 0, . . . } = {0, . . . , 0, e - iHNt <p m , 0, . . . } 

where Hn is the mean field Hamiltonian defined in ([5]). What did we gain then by switching 
to a Fock space representation of the many boson system? The gain is that now we have 
much more freedom in the choice of the initial state. We will make use of this freedom by 
considering so called coherent states as initial data. 

For / G L 2 (M 3 ), we define the Weyl operator 

W(/) = exp(o*(/)-a(/)) . 
The coherent state with wave function / is then defined as 

W (/)n = e HI/IIV2^(^^ n = e HI/llV2( lj/> /^ > . ..) . (9) 

n>0 L V2! J 

Since W(f) is a unitary operator on J 7 , coherent states are always normalized. It is clear 
from (jSJ) that coherent states do not have a fixed number of particles; instead they are given 
by a linear combination of states with all possible number of particles. Still, one can ask 
what is the expected number of particles in the coherent state W(/)f2. A simple computation 
shows that 

(w(m,Afw(f)n) = |M| 2 . 

This follows from the fact that Weyl operators act on creation and annihilation operators 
by simple shifts, in the sense that 

W*(f)a*(g)W(f) = a*(g) + (f,g) 

W*(f)a(g)W(f)=a(g) + (g,f). [ j 

The property ()10p also implies that coherent states are eigenvectors of all annihilation oper- 
ators, i.e.: 

a(g)W(m = (g,f)W(f)n 

for every f,g € L 2 (M 3 ). This important observation is the main reason for the nice algebraic 
properties of coherent states which, as we will see below, significantly simplify the analysis 
of their time evolution. 

We study next the evolution of coherent states. In order to recover the mean field regime 
discussed in the introduction, the number of particles must be related with the parameter N 
appearing in the Hamiltonian (JS} (equivalent to ©)• Since we want the initial state to be 
coherent, we cannot fix the number of particles to be equal to N. Still, we can require N to 
be the expected number of particles. Hence, we fix ip £ L 2 (IR 3 ) with \\(f\\2 = 1, and consider 
the initial coherent state 

y N = w(y/Ntp)n. 

We study its time evolution 

y Nit = e- iHNt ^ N = e- iHNt W(^Ncp)Q 



6 



as generated by (jSJ). In particular, we would like to compute the reduced densities T Nt 
associated with ^>N,t- Let us consider for example the one-particle reduced density. It turns 

out that the kernel T^ t (x,y) of the one-particle reduced density can be expressed as 

T$ t (x,y) = — (^N,t,ala x ^ N}t ) 

1 — W 

= — (fr iHNt W{VN<p)tt, a* y a x e~ iHNt W(VN(p)nj . 

Of course, for t ^ 0, ^Nt is n °t a coherent state. Nevertheless, because of the mean- field 
character of the interaction, in the limit of large N we can expect it to be approximately 
coherent, with an evolved one-particle wave function y/N<pt, where (ft solves the Hartree 
equation @. In other words, we may expect that 

y Nt = e - iHNt w(VNip)n ~ w(VN(p t )n . 

Recalling that coherent states are eigenvectors of the annihilation operators, it seems appro- 
priate to expand the operator valued distributions a x ,a* on the right hand side of (fTT|) around 
their mean field values yNiptix) and \HSlTp t (y). Writing a x = yNipt(x) + (a x — yNcp t (x)) 
and similarly for a*, and inserting in we find 

= 1 (n, W*(VN<p)e m ^ (a* y - VNlp t (y))(a x - Vn^x)) e -^^(^)fi) 




Observe here that (pt{x)lp t {y) is the kernel of the orthogonal projection \ipt){tpt\ onto the 
solution of the Hartree equation @. To establish the convergence of the full many-body 
evolution towards the limiting Hartree dynamics, we need therefore to bound the error terms 
on the r.h.s. of (fT2|h 

Next, we make use of an idea introduced in the related context of the classical limit 
of quantum mechanics in [5]. We observe namely that, from ()10p . the fluctuations of the 
operator valued distributions a*,a x around their mean field values can be expressed as 

(a* y - VNWM) = W{VNift) a* y W*(VN Vt ) 
(a x - y/N<pt(y)) = W{VN^t) a x W*(y/N<p t ) • 

From (|12p . we conclude that 




7 



Introducing the fluctuation dynamics described by the two-parameter group of unitary trans- 
formations 



U(t; s) = W(VNip t )e- iHN{t - s) W*(VNip s ), for all t,sGl, (13) 
we obtain the identity 

= ^(nM*(t;0)a* y a x K(t;0)n) (14) 



Let us focus on the first error term on the r.h.s. of (|14p . Since, by ([7]), creation and 
annihilation operators are bounded with respect to the square root of the number of particles 
operator, the smallness of this term follows if we can show uniform (in N) estimates for the 
growth of the expectation of M with respect to the dynamics U(t; 0). To obtain such bounds, 
we observe that the fluctuation dynamics satisfies 

id t U(t; s) = C N (t)U{t; s) with U(s; s) = 1 for all s E M 

with the time-dependent generator 

£jv(i) = j dxV x a* x V x a x + j dx (V * \ip t \ 2 ){x)a* x a x 

+ J dxdy V(x - y)ip t {x)Tp t (y) a* x a y 

+ dxdyV(x-y)(^ t (x)ipt(y)a x a* y +Tp t (x)Tp t (y)a x a y ) (15) 



1 



+ -7= j dxdyV(x - y) a* (<p t (y)a y + ft(y)a* y ) a x 



^ J dxdy V(x - y)a* x a* y a y a x . 



+ N 

In contrast with Hn, the generator C]\[(t) (in particular, the terms on the third and fourth 
line) does not commute with Af. As a consequence, the fluctuation dynamics U (t; s) does 
not preserve the number of particles. This is hardly surprising, since fluctuations around the 
mean field solution have to grow during the time evolution. Although the expectation of the 
number of particles is not preserved under U{t\ s), it is possible to show uniform bounds for 
the growth of M and any of its power. The following result was proven in [7]. 

Proposition 2.1. Suppose that there exists a constant D > such that the operator inequal- 
ity V 2 {x) < D{1 — A) holds true. Suppose moreover thatU(t;s) is defined as in <f73|). Then, 
for every k S N there exist constants C, K > with 



>, U*(t; s) (AT + l) k U{t- s)ip) < Ce^-^V, {N + l) 2k+2 V) 
for all t, s E M. 

Remark. The operator inequality V 2 (x) < D(l — A), meaning that 

j dxV{x) |v?(x)| 2 <D J dx [|V(^(x)| 2 + \ip(x)\ 2 ] = D\\<p\\ 2 H1 , 



S 



for all (p G L 2 (M 3 ), is satisfied, because of Hardy's inequality, for potentials with Coulomb 
type singularities V(x) = ±l/|x|. 

Proposition 12.11 immediately implies that the first error term on the r.h.s. of (|14|) is of 
the order 1/N, for any fixed t G R. With some more work one can show the same estimate 
also for the last two terms on the r.h.s. of (|14p . As a consequence, one obtains convergence 
towards the Hartree dynamics. The details of the proof of the next theorem can be found in 

Theorem 2.2. Suppose that there exists a constant D > such that the operator inequality 
V 2 (x) < D(l — A) holds true. Suppose that, at time t = 0, = W(y/~Nip)Q, for some 
<p G ^(R 3 ). Let V Nt = e~ mNt ^ N , with the Hamilton operator and let F$ t be the 



one-particle density associated with ^N.t- Then there exist constants C, K > with 



tr 



t n\ ~ \<Pt)(<Pt 



< 

- N 



where (ft is the solution of the Hartree equation with initial data y>t=o = V 3 - Similar 
bounds can also be proven for the k-particle reduced density, for any fixed k G N. 

Remark. The same result can be obtained, with exactly the same techniques, for initial 
states of the form ipx = W(y/N(p)ip, for a ip G T with (i^,Af 2 ip) < 1 (independent of N). 

This approach to the study of the evolution of coherent states not only implies the 
convergence towards the mean field Hartree dynamics. Instead, it also establishes the form 
of the fluctuation dynamics in the limit of large N. Formally, the generator C^it) converges, 
as iV — > oo, towards the Fock-space operator 

j &xV x a* x V x a x + J dx (V * \(p t \ 2 ){x)a* x a x 

+ J dxdy V(x - y)ip t {x)Tp t (y) a* x a y 

+ / dxdyV(x -y)((pt{x)(p t {y)alal+Tp t (x)Tp t {y)a x ay) . 



One can expect, therefore, that the fluctuation dynamics U(t;s) converges, as iV — > oo, 
towards the limiting dynamics Uoo{t; s), defined by 

idtUooit; s) = £ 00 (t)^/ 00 (t; s) with U^s; s) = 1 for all s G R . 

Since £oo(^) is a time-dependent unbounded operator, the definition of the limiting dynamics 
Uocit; s) generated by £oo{t) is not at all trivial. The existence of Uoo{t; s), and the (strong) 
convergence U(t; s) — > Uoo(t; s) were rigorously established in [4J, making use of appropriate 
approximations of Coc (t) . 

Since the generator £oo(t) is a quadratic expression in creation and annihilation operators, 
it turns out that the limiting dynamics U^t; s) can be described as a so called Bogoliubov 
transformation. For /, g G L 2 (R 3 ), we define on T the linear combination of creation and 
annihilation operators 

A(f,g) = a*(f) + a(g). 
By definition, A is linear in both / and g. Observe that 



(A(f,g)Y =A(Jg,Jf)=A 



J\(f 
J M g 
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where J is the antiunitary operator on L 2 (R 3 ) defined by Jf = f. A simple computation 
shows that, in terms of the operators A(f, g), the canonical commutation relations © assume 
the form 

fx W i o w / 2 

ffl J ' V -1 J V 52 

where (., .) denotes the standard inner product on L 2 (M 3 ) © L 2 (R 3 ). A Bogoliubov trans- 
formation is a linear map 9 : L 2 (R 3 ) © L 2 (R 3 ) -»■ L 2 (R 3 ) © L 2 (M 3 ) with the properties 



Wi,3iM(/2,32) 



and 






J 


J 




' 1 










J 
J 



1 

0-17" V o -l 



(16) 



(17) 



Eq. (fT6|) guarantees the preservation of the relation between A and its adjoint. Eq. (fTT|) . 
on the other hand, guarantees the preservation of the canonical commutation relations. 
It is simple to check that every Bogoliubov transformation can be expressed through the 
operator-valued matrix 

■ U JVJ\ 
V JUJ ) (18) 

where U, V : L 2 (R 3 ) -»• L 2 (R 3 ) are s.t. U*U - V*V = 1 and U*V - V*U = (notice that 9 
is not unitary, unless V = 0). Every Bogoliubov transformation defines therefore a new set 
of creation and annihilation operators. It is interesting to ask when the new representation 
of the canonical commutation relation is unitary equivalent to the one given by the original 
operators a*(f), a(f). It turns out that this is the case if and only if the operator V appearing 
in (|18|) is Hilbert-Schmidt (Shale-Stinespring condition). 

The statement that the limiting fluctuation dynamics Uoo(t;s) acts as a Bogoliubov 
transformation has to be understood as follows. There exists a two parameter family of 
Bogoliubov transformation 0(t; s) : L 2 (M 3 ) © L 2 (M 3 ) ->■ L 2 (R 3 ) © L 2 (M 3 ) such that 

U* ca (t;s)A(f,g)U 00 (t;s) = A(e(t; S )(f,g)) (19) 

for every f,g £ L 2 (R 3 ). Formally, 8(t;s) is given by the solution of 

A -Bt 



id t 9(t;s)=[ B ^ )9(t;s) 



with initial condition 9(s; s) = 1, and with D t , B t : L 2 (R 3 ) — > L 2 (R 3 ) defined by 

D t f = -A/ + (V * \<p t \ 2 )f + (V * lp t f)<pt 
Btf = (V*Tp t f)Tp t . 

We will not make use of this formal characterization of 9(t; s). The important observation 
is that the time evolution Uoc(t;s), which is in principle a two-parameter family of unitary 
transformation on the large Hilbert space J- can be completely described in terms of the 
family 9(t; s) operating on the much smaller space L 2 (M 3 ) © L 2 ( 



So far we considered the evolution of initial coherent states. Next, we turn our attention 
back to initially factorized (or approximately factorized) iV-particle states, as those consid- 
ered in the introduction. To this end, we notice that, for any (p £ L 2 (R 3 ), we can write the 
factorized state 

{o, . . . , o, <p m ,o, ...} = d N p N w(VNtp)n 
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where P/v denotes the orthogonal projection onto the sector with exactly N particles, and 
where d n is a normalization constant; a simple computation shows that djy ~ iV 1 / 4 . The 
reduced one-particle density associated with the evolution of the factorized initial data is 
given by 

/ e _ ffljyt (°*(g))" n> Q * Qye - ffl ^ t P jV W(v / iVy)0 N ) 



iV \ VM x y v y / 



because P/v commutes with and with a*a x , and because Pjy (a* (y)) ^ = (a*(^)) 7V ri. 
Letting 

e = div^*(Vjvy) (Q ^ )Ar o 

we find 

7#J(*,v) = ^(e,W(t;0)(a* + x/jV? t (x))(aj, + VN<p t (y))U*(t;0)n 

= Tp t (x)<pt(y) + 1 0)a* x a y U(t; 0)Q) (20) 

+ ^ (£, W*(t; 0)a,W(t; 0)17) + ^ )<W(t; 0)17) 



where U(t;s) denotes the fluctuation dynamics defined in (|13|) . From this formula, we 
conclude that, similarly as for initial coherent states, proving the convergence towards the 
Hartree dynamics reduces to the problem of obtaining uniform bounds for the growth of the 
product 

\(t,U*(t;0)ArU(t;Q)n) \ . (21) 

The only difference compared to the case of coherent initial data is that now, on the l.h.s. 
of the product, we have the vector £ instead of the vacuum. Since ||£|| = d N ~ iV 1 / 4 , it 
seems now more difficult to get estimates uniformly in N. It turns out, however, that when 
restricted to sectors with small number of particles, £ is an order one vector, in the sense 
that 

IKAA + irtH <1 

uniformly in N. For this reason, 

\(Z,U*(t;0)MU(t;0)n)\ < IKAf+l)- 1 ^! \\(M+l)K*(t;0)MU(t;0)n\\ 

<e A 'l*l || {M+lfU{t; 0)n|| (22) 
< e 2*|t||| (A f +1) 7 fi || < e 2*|t| 

applying twice Proposition [2TTJ As a corollary of (|22p . with some additional work needed to 
bound the last two terms on the r.h.s. of (|20p . we obtain the convergence of the evolution of 
initially factorized data towards the Hartree dynamics, with an explicit bound on the rate 
of convergence. The proof of the following theorem was found in [2j, optimizing ideas from 
[7j (in [7J, the error for factorized initial data was shown to be at most of the order iV -1 / 2 , 
for every fixed t E R). 
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Theorem 2.3. Suppose that there exists a constant D > such that the operator inequality 
V 2 (x) < D(l — A) holds true. Let ipN.t = e~ lHNt <p® N with the Hamilton operator and 
for some ip E H 1 (M 3 ). Let 7^ be the one-particle reduced density associated with ipN,t- 
Then there exist constants C, K > with 



tr 



\<Pt)(<pt 



< 



N 



Similar bounds holds for the k-particle reduced densities as well. 

Remark. The same result can be obtained for initial data of the form ip^ = <InPnW( yN tp)tp, 
for an arbitrary i/j E T with {^,M m tp) ^ 1 for some m E N large enough, and where 
d/v — A rl / 4 is chosen s.t. ipN is normalized. 



3 A probabilistic setting 

In this last section, we formulate the convergence towards the mean field Hartree dynamics 
in a language more common in probability theory. To this end, we consider an TV-particle 
system, described by a permutation symmetric wave function ip]\f E L 2 (M? N ). We consider, 
moreover, a self adjoint operator O acting on L 2 (M. 3 ). For j = 1, ... ,JV, we define O^' to 
be the self-adjoint operator acting as O on the j-th particle and as the identity on the other 
(N — 1) particles. Every O^' can be thought of as a random variable assuming different 
values with different probabilities. Through the spectral theorem, the wave function ipjy 
determines the law of the random variables O^'. 

At time t = 0, we assume that the system is described by a factorized wave function 
^tv = <P® N 1 for some ip E L 2 (R 3 ). The operators j = 1, . . . , N, define then a sequence 
of independent and identically distributed random variables, with a common distribution 
determined by ip. We obtain immediately the (weak) law of large numbers, stating that, for 
any 5 > 0, 

1 N 

3=1 



> 5 







as N — > 00. We also have a central limit theorem, stating that, in distribution, 



1 N 

-J2(oW-(<p,O<p))^N(0,o-> 
3=1 



where N(0, a 2 ) denotes a centered Gaussian random variable, with the variance 

o 2 = (<p,0*<p)-(<p,0<p) 2 . 



Now, let us consider the evolution of the many body quantum system, as generated by 
the mean field Hamiltonian (|2|) . The wave function describing the evolved system is given by 
ipN,t = e~ lHNt p® N . For t ^ 0, ipN,t is n °t factorized; hence, the random variables O^' are 
not independent. Still, the results presented above imply that correlations are small in the 
limit of large N, at least in the sense of the reduced densities. It seems therefore natural to 
ask whether law of large numbers and central limit theorem are still valid at t 7^ 0. It turns 



12 



out that the convergence of the reduced density, as stated in Theorem 12. 3\ easily implies the 
(weak) law of large numbers. In fact, letting O = O — ((ft, Oipt), we find 



N 



N ^ 



> S < 



1 



' N 
a=l 

^tr|^)(^| 2 (d®0) = 




because, by definition of O, tr \ip t }(y>t\0 = ((ft, Otpt) = 0. 

What about the central limit theorem at t 7^ 0? The answer to this question is given in 
the next theorem, which was recently proven in [TJ. 



Theorem 3.1. Suppose that V 2 (x) < D(l — A) for a constant D > 0. Let i/^N,t = e 



with Hjy given in H|) and for ip € H l (E?). Let O be a bounded self-adjoint operator on 
L 2 (M 3 ), with ||VO(l — A) _1 / 2 || < co. Then, for every t € R, with respect to the evolved wave 
function ipN,t, the random variable 



1 N 

i=i 



converges in distribution, as N — )• oo, to a centered Gaussian random variable with variance 
(0(i;O) (0|ft,0^) ,e(t;0) (0^,0^) > 



e(t;0) {0<pt,0<p t ),-j= (<p,lp) 

where 9(t;0) is the time- dependent Bogoliubov transformation introduced in Fffl\). 

Hence, for t ^ 0, the central limit theorem is still valid but the variance of the limiting 
Gaussian variable is changed. While on the level of the law of large numbers there is no 
difference between ipN,t and the product ipf N , this is not true for the central limit theorem; 
although in both cases the fluctuations are gaussians, their variance is not the same. 
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